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We argue that the motion of vacancies in a pinned vortex lattice may dominate the contri- 
bution to the Hall effect in an appropriate parameter regime for a superconductor. Based 
on this consideration a model is constructed to explain the anomalous Hall effect with- 
out any modification of the basic vortex dynamic equation. Quantitative predictions are 
obtained. Present model can be directly tested by an observation of the vacancy motion. 
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It is generally believed that the motion of vortices in 
a superconductor is responsible for the apparence of re- 
sistance. Indeed, there is a well and successful account 
of dissipative phenomena in superconductors based on 
this ideaQ. However, the ubiquitous occurrence of the 
anomalous Hall effect in the mixed state of both con- 
ventional and oxide superconductors, the sign change 
of the Hall resistance below the superconducting tran- 
sition temperature, seems to be inexplainable by stan- 
dard models for vortex dynamics in superconductors^) . 
The purpose of the present paper is to show that it is 
nevertheless possible to understand the anomalous Hall 
effect without either a modification of vortex dynamics 
or assuming two types of carriers(or more complicated 
electronic band structures). The essential idea in the 
present model is the domination of the vacancy motion 
in a pinned vortex lattice under a suitable condition. 
Contrast to previous ones in this model both the strong 
interaction between vortices, which leads to the many- 
body correlation effect, and the strong pinning effect 
are crucial. In the following we present the arguments 
leading to the model, and discuss its predictions and 
experimental verifications. For simplicity, we will con- 
sider an isotropic s-pairing superconductor with one 
type of charge carriers in two dimension. In this situa- 
tion vortices (or straight vortex lines) can be viewed as 
point particles. 
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We begin with the discussion of vortex dynamics. 
The equation of motion for a vortex generally takes 
the form of the Langevin equation similar to that of a 
charged particle in the presence of a magnetic field: 

m v£ = qv-^hd (v s - r) x z - rjr + F pin + f , (1) 

with an effective mass m v , a pinning force F P i„, a vor- 
tex viscosity 77, and a fluctuating force f. The vis- 
cosity is related to the fluctuating force by the usual 
fluctuation-dissipation theorem. In Eq.(l) q v = ±1 de- 
scribes the vorticity, h is the Planck constant, p s the 
superfluid electron number density at temperature T, 
d the thickness of the superconductor film, and z the 
unit vector in z-direction. Normal current is assumed 
to be zero and a uniform temperature in the whole film 
is also assumed in Eq.(l). The term associated with 
the superconducting electron velocity v s and the vor- 
tex velocity r at the right side of Eq.(I) is the so-called 
Magnus force. It has been shown that under the gen- 
eral properties of a superconductor it must take such 
a form||. Closely related but different phenomenolog- 
ical equations for vortices have been discussed in the 
literature Eq.(l) is the starting point in recent stud- 
ies of vortex tunneling in superconductors ||. 

We first consider the measured electric field gener- 
ated by a steady motion of vortices and the resulting 
Hall effect. In a steady state, if there is no pinning 
force, i.e. F p i„ = 0, from Eq.(l) the velocity of the 
vortices is 
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According to the Josephson relation, the motion of vor- 
tices produces phase slippages, therefore generates an 
electrochemical potential difference in a superconduc- 
tor, which is the reaction force to the superconducting 
electrons due to the vortex velocity part of the Magnus 
force in Eq.(l). This potential difference corresponds 
to the measured electric field E. For a vortex density 
n with an average velocity r = v; , we have 



E = -q v — tivi x z 
2e 



(3) 



If the density of moving vortices is equal to the vortex 
density induced by the applied perpendicular magnetic 
field B = Bz, n = \B\/$ with <J> = hc/2\e\ the flux 
quantum of Cooper pairs, we have the usual expression 
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Therefore, in the steady state and in the absence of 
pinning, all vortices move with the same velocity and 
we have the longitudinal resistivity p xx as 
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and the transversal(Hall) resistivity p yx as 
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Here the current density J = ep s v s is along x-direction. 
The relation between vorticity q v , the charge sign of 
carriers, and the direction of the magnetic field is 
q v e/\e\ — B/\B\. It is clear that from Eq.(6) the Hall 
effect in the superconducting state has the same sign 
as that of the normal state. This, and also those sim- 
ilar conclusions based on the classical models for the 
vortex dynamics in superconductors, are in apparent 
contradiction with many experimental observations g, 
the so-called anomalous Hall effect. 

There is a serious problem related to the many-body 
correlation between vortices and the pinning effect in 
the using of the above simple minded result, Eq.(6), for 



real superconducting materials. The many-body cor- 
relation due to the strong and long range interaction 
between vortices gives arise to the Abrikosov lattice in 
the mixed state of a superconductor. However, even 
after the consideration of the strong many-body effect, 
if there were no pinning for vortices, the whole vor- 
tex lattice would move together under the influence of 
an externally applied current in the same manner as 
that of independent vortices. Hence one would get the 
same sign of the Hall effect in both superconducting 
and normal states. This will be changed in the pres- 
ence of pinnings. If we introduce some strong pinning 
centers into the film and consider the edge pinning, the 
vortex lattice will be pinned down. In such a situation 
the motion of the vortex lattice is made possible by var- 
ious kinds of thermal fluctuations. We will argue below 
that at low temperatures the dominant contribution to 
the motion is due to vacancies. We will advance our 
argument in the following way. First we will show that 
the creations and motions of vacancies and interstitials 
are energetically favored over vortices hopping into and 
out of pinning centers as well as other motions. Then 
we will argue that vacancies are even more favourable 
than interstitials. 

For two vortices separated by a distance r, which is 
less than the effective magnetic screening length Aj_ = 
\ 2 L /d(d < Xl, X± = Xl if d > Xl) but greater than £ , 
the interaction potential ispfl 
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Here X\ — m*c 2 /8Trp s e 2 is the London penetration 
depth, m* the effective mass of a Cooper pair, and 
£o the coherence length of the superconductor. The 
energy scale eo sets both the scale for the strength of 
vortex interaction and the scale for the strength of a 
strong pinning center. The interaction potential be- 
tween a dislocation pair separated by a distance larger 
than the lattice constant ao is given by[Q 



V d (r) = _^ eo ln- = 2eln- 
2V37T ao ao 



(8) 



The energy scale e for the dislocation interaction here 
is about 20 times smaller than eo for the vortex inter- 
action and pinning centers. It is therefore energetically 
much more favorable to have dislocation pairs in the 
lattice. Hence for temperature T << eo we can ignore 
the contribution from the vortices hopping out of pin- 
ning centers (and also the creation of vortex- antivortex 
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pairs). The vortex lattice is then effectively pinned at 
such temperatures. The precise calculation of the va- 
cancy e v and interstitial energy is difficult, because 
the detailed information on the scale of lattice constant 
is needed. However, because vacancies and interstitials 
can be viewed as the smallest dislocation pairs ||, we 
immediately have the estimated energy scale for e„ as, 
by putting r ~ 2ao in Eq.(8), 

e v ~ 2e = — ^- ( -^-) d. (9) 
2V37T V 47r -W 

There might be a weak dependence on the magnetic 
field through the vortex lattice constant, which is ne- 
glected here. It is clear from the above analysis that 
vacancies and interstitials have the lowest excitation 
energy scale. 

Now let us compare the creation energy of a vacancy 
and that of an interstitial. At low magnetic fields the 
vacancy formation energy is lower than that of inter- 
stitials, because of the strong short range and repul- 
sive nature of the interaction between vortices. This 
has been shown theoretically |9j. Therefore vacancies 
will dominate over interstitials at low temperatures and 
low fields, which has been observed experimentally . 
Similar phenomenon has also been observed in other 
crystalline structures At high magnetic fields the 
vortex interaction becomes of long range comparing 
with the vortex lattice constant. The short range be- 
havior which determines the formation energy differ- 
ence between vacancies and interstitials becomes rela- 
tively unimportant. Although one expects a zero for- 
mation energy difference for an infinite range interac- 
tion, vacancies may still have a lower formation energy 
for a finite range interaction. This is based on the plau- 
sible expectation that the formation energy difference 
between vacancies and interstitials increases monoton- 
ically to zero with the interaction range. We will not 
be surprised if the anomalous Hall effect is more pro- 
nounced at lower magnetic fields. 

We examine its experimental consequencies if the va- 
cancy motion dominates. Effectively, the motion of a 
vacancy in the pinned vortex lattice behaves as an an- 
tivortex, that is, with a vorticity —q v under the action 
of an applied supercurrent. This leads us to our main 
conclusion that at low enough temperatures the sign 
of the Hall resistance is different from its sign in the 
normal state. Quantitatively, the motion of vacancies 
is governed by an equation of the same form as Eq.(I), 
and can be viewed as independent particles moving in 
the periodic potential formed by the vortex lattice and 
a random potential due to the residue effect of pin- 



nings. The potential height of the periodic potential 
as well as that of the random potential is presumably 
the order of e v . Assuming the vacancy density n v in 
a steady state, following the same procedure as from 
Eq.(I) to Eqs.(5-6), the longitudinal resistivity is then 

h r] ef f p s dh/2 n v 
Pxx ~ 2e* t fS + (Psdh/2y p s ' UUj 

and the Hall resistivity 

h (p s dh/2) 2 n v 
Pyx ~ 9v 2e> if eff + (p a dh/2y p s ■ {U) 

It is evident that the sign for p yx in Eq.(II) is opposite 
to that in Eq.(6). Here ry e // is the effective vacancy vis- 
cosity, related to the vacancy diffusion constant in the 
periodic potential due to the vortex lattice by the Ein- 
stein relation between the diffusion constant and the 
mobility. We expect that r] e ff = tjq e a<L v/K B T w j£ n a 
a numerical factor of order unity at low temperatures. 
In the low temperature limit, superfluid electron num- 
ber density p s approaches a constant, vacancy density 
n v — uq e - be v/kisT ma y b e exponentially small(6 = 1 
for the thermally activated vacancies, and 6 = for the 
pinning center induced vacancies), and the effective va- 
cancy viscosity rj e ff can be exponentially large. In this 
limit both longitudinal and Hall resistances vanish ex- 
ponentially, and we obtain a scaling relation between 
the Hall and longitudinal resistivities as 

Pyx =Apx X , (12) 

with A = -q v (p s dh/27] ) b /( a+b )(2e 2 p s /hn ) a /( a+b ) and 
the power v = (2a + b)/(a + b) varies between I and 
2 depending on the detail of a sample which deter- 
mines the numerical factors a and b. It is inter- 
esting to note that the present range of v covers 
the values obtained previously by completely different 
approaches Jll|. We also note that if the vortex (va- 
cancy) tunneling process^) is allowed, there will be a 
finite Hall resistance at zero temperature. 

The linear film thickness d dependence of the va- 
cancy formation energy e v , Eq.(9), can also be tested 
experimentally. We expect this linear dependence to 
hold for thin enough films. For a thick film a vortex 
line will bend in the z-direction, and the straight line 
assumption in reaching Eq.(9) is no longer valid. This 
suggests a critical thickness d c such that Eq. (9) is valid 
for d < d c and e v will be saturated for d > d c . This 
critical thickness may reflect the many-body correla- 
tion along the direction of the magnetic field and is 
likely a sensitive function of both the magnetic field 
and the temperature. 
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Now we discuss the qualitative features of the model. 
In the above picture, we need the vortex lattice to de- 
fine vacancies, and a sufficiently strong pinnings to pre- 
vent the slidding of vortex lattice in order to obtain a 
maximum contribution of vacancies. However, there is 
no need for a whole lattice structure. Sufficiently large 
local structures, like lattice domains, will be enough to 
define vacancies. Therefore one may even have vacancy 
like excitations in a vortex liquid state, where large lo- 
cal orderings exist. Whether or not this is also true 
for a vortex glass state depending on its randomness. 
On the other hand, if the pinning is too strong, for ex- 
ample, the pinning center density is much larger than 
the vortex density, vortices will be individually pinned 
down and the local lattice structure required for the 
formations of vacancies and interstitials will be lost. 
Therefore we only expect to see the anomalous Hall ef- 
fect in a suitable range of pinnings and magnetic fields, 
that is, for Bi < \B\ < B u with the lower and upper 
critical fields determined by pinnings. 

As shown by Eqs.(7-9) a large London penetration 
depth reduces down all the relevant energy scales, and 
energetically it is more favorable to observe the anoma- 
lous Hall effect. This may explain the experimental 
observation that it is difficult to observe the anoma- 
lous Hall effect in a clean conventional superconduc- 
tor, but easy in a dirty superconductor as well as in 
a high temperature superconductor where the London 
penetration depth is large 0. At this point, we wish to 
point out the important difference between the pinning 
centers for vortices and the mean free path of electrons: 
The vortex dynamics is sensitive to pinning, not to the 
mean free path. There is no obvious relationship be- 
tween them. 

If the temperature is increased towards the supercon- 
ducting transition temperature, all energy scales be- 
come small, and so does the difference between e v and 
€i. Other contributions, such as those of quasiparticles, 
the slidding of vortex lattice, and vortices hopping out 
of pinning centers, become dominant. Those contribu- 
tions have a different sign for the Hall resistance, but 
the same sign for the longitudinal resistance as vacan- 
cies. Therefore one expects a sign change in the Hall 
resistance below the transition temperature. 

It is important to check whether or not the present 
model can also explain experiments concerning vortex 
motion other than the Hall and longitudinal resistivity 
experiments. We will mention here the Nernst effect. 
Under the driving of a temperature gradient the force 
felt by a vacancy is opposite of the force felt by an 
interstitial or a vortex in direction but equal in magni- 



tude. Then the Nernst effect due to vacancies has the 
same sign as that of vortices or interstitials. There- 
fore our model gives that in the anomalous Hall effect 
regime there is no sign change for the Nernst effect, 
and furthermore, the Nernst effect is more pronounced 
because of the large contribution due to both vacan- 
cies and interstitials. This is in agreement with the 
experimental observations . 

In conclusion, a model based on the motion of va- 
cancies in a pinned lattice is proposed to explain the 
anomalous Hall effect. No modification of the basic 
vortex dynamics equation is needed. Quantitatively it 
leads to an exponential tail and the scaling relation at 
low temperatures, and no sign change for the Nernst 
effect. For thin enough films the activation energy in 
the low temperature limit has a linear film thickness 
dependence. Present model provides a framework to 
understand relevant experiments, and further experi- 
ments are needed to test it. The most direct way might 
be the observation of the vacancy motion. 
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